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We investigate the possibility of realizing quantum anomalous Hall effect in graphene. We show 
that a bulk energy gap can be opened in the presence of both Rashba spin-orbit coupling and an 
exchange field. We calculate the Berry curvature distribution and find a non-zero Chern number 
for the valence bands and demonstrate the existence of gapless edge states. Inspired by this finding, 
we also study, by first principles method, a concrete example of graphene with Fe atoms adsorbed 
on top, obtaining the same result. 
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Graphene is a fascinating material not only because 
of its extraordinary electrical, magnetic, and mechanical 
properties, but also because of its Dirac-cone structure in 
its low energy spectrum, which involves three binary in- 
ternal degrees of freedom: spin, sublattice and valley jll- 
y] . Manipulating these degrees of freedom may turn out 
to be the most rewarding means of controlling the prop- 
erties of graphene and related materials. For example, by 
breaking the AB sublattice symmetry, an energy gap can 
be opened at the Dirac points [5-8], leading to the valley- 
contrasting Hall, magnetic and optical responses [9]. The 
intrinsic spin-orbit coupling (SOC) is proposed to give 
rise to the quantum spin Hall effect [10, 11], a quantized 
response of a transverse spin current to an electric field. 

However, the intrinsic SOC in graphene is proved to be 
too weak to realize the quantum spin Hall effect under 
present experimental conditions 12|, ll3|- On the other 
hand, the extrinsic Rashba SOC from breaking the mir- 
ror symmetry of the graphene plane tends to destroy this 
effect [10]. Rashba SOC can be very large for graphene 
grown on a substrate [ij, ll5|- For example, a 225 meV 
of Rashba spin splitting has been observed for graphene 
grown on Ni [15]. We expect that Rashba SOC should 
also be present for graphene on insulating substrate, or 
with dilute adsorbates. In the following, we show that al- 
though Rashba SOC is detrimental to the quantum spin 
Hall effect, it helps to realize another important topolog- 
ical phenomenon: the quantum anomalous Hall effect. 

This effect is signaled by a quantized charge Hall con- 
ductance for an insulating state. Unlike the quantum 
spin Hall effect, the quantum anomalous Hall effect re- 
quires the breaking of time reversal symmetry. Unlike 
the quantum Hall effect, which is induced by a strong 
magnetic field, the quantum anomalous Hall effect relies 
on the internal magnetization and SOC. The quantized 
value of the Hall conductance is related to a bulk topo- 
logical number and is robust against disorder and other 
perturbations [16|, ll7| . Although there have been many 



theoretical studies of this unusual effect J18l-l22| , no ex- 
perimental observation has been reported so far. 

In this Letter, we predict that the quantum anoma- 
lous Hall effect can be realized in graphene by introduc- 
ing Rashba SOC and an exchange field. We also con- 
sider a concrete example with transition metal atoms 
{e.g., Fe) adsorbed on top of graphene. The resulting 
broken structural symmetry gives rise to a Rashba type 
SOC whereas the hybridization between the carbon tt 
state and the 3(i-shell states of magnetic atoms produces 
a macroscopic exchange field. The Hall conductivity is 
found to be dyx = 2e^ /h when the Fermi level lies in 
the bulk gap. From first principles calculations, we find 
that a bulk gap ~ 5.5 meV can be opened by adsorbing 
Fe atoms on top of graphene, which is readily accessible 
under current experimental conditions. 

In the tight-binding approximation, the Hamiltonian 
for graphene in the presence of Rashba SOC and ex- 
change field can be written as: 

{ij)a {ij)a(3 

+ ^X^4c.^^ci« (1) 

ia 

where q(Q) is the electron creation (annihilation) op- 
erator on site z, and a are the Pauli matrices. The 
angular bracket in (i,j) stands for nearest neighboring 
sites. The first term is the usual nearest neighbor hop- 
ping term. The second term is the Rashba SOC with 
coupling strength t^o, and dij represents a unit vector 
pointing from site j to site i. The third term corresponds 
to a uniform exchange field [23]. Throughout this Letter, 
we measure the Fermi level, Rashba SOC, and exchange 
energy in units of the hopping parameter t. 

By transforming the Hamiltonian of Eq.([T]) into a 4 x 4 
matrix H{k) for each crystal momentum, the band struc- 
ture of bulk graphene can be numerically obtained by 
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FIG. 1: Evolution of band structures of bulk graphene along 
the profile of ky = 0. Arrows represent the spin directions, (a) 
Pristine graphene: spin-up and spin-down states are degen- 
erate; (b) When only exchange field A = 0.4 is applied, the 
spin-up/spin-down bands are upward /downward lifted with 
the four bands crossing near K and K' points; (c) When only 
Rashba SOC tso = 0.1 is present, the spin- up and spin-down 
states are mixed around the band crossing points; (d) When 
both exchange field A = 0.4 and Rashba SOC tso = 0.1 are 
present, a bulk gap is opened and all four bands become non- 
degenerate. 



diagonalizing H{k). Fig(T] shows the evolution of the 
band structures with Rashba SOC tgo and exchange field 
A. Panel (a) plots the band structure of the pristine 
graphene with Dirac cones centered at K and K^ points 
in the reciprocal space. Due to spin degeneracy, K and 
K' points are four- fold degenerate, while other points 
are doubly-degenerate. When only the exchange field is 
applied, the spin-up (spin-down) bands are pushed up- 
ward (downward) as shown in panel (b). When only 
Rashba SOC is turned on, spin-up and spin-down states 
are mixed around the band crossing points and spin de- 
generacy is lifted. A bulk gap is opened when both 
Rashba SOC and exchange field are present, and the four 
bands become completely non-degenerate. The presence 
of such a bulk gap indicates an insulating state. As we 
argue below, this insulating state is topologically non- 
trivial with gapless chiral edge states and exhibiting a 
quantized charge Hall conductance. 

The gapless chiral edge states can be clearly seen from 
the energy spectrum of graphene ribbons. Fig 121 (a) plots 
the band structure of graphene ribbons, with zigzag edges 
and 800 atoms across the width. The parameters used 
here are tgo = 0.2 and A = 0.18. One can easily dis- 
tinguish the gapless edge states from the bulk states. 
For a given Fermi level in the gap, there exists four dif- 
ferent edge states labeled as A, 5, C, and D. From 
v{k) = ^ QJ^\ one can find that states A and B {C and 
D) propagate along the same —x {-\-x) direction. Panel 
(b) plots the wave function distributions l^/^p of the four 
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carbon atom sequence across the ribbon 

FIG. 2: (Color online) (a) Energy spectrum of zigzag-edged 
graphene ribbons with tso = 0.2 and A = 0.18. The Fermi 
level E = 0.06 corresponds to four different edge states A, 
B, C, and D. a is the lattice constant, (b) Wave function 
distributions l^^l^ across the width for the four edge states 
(only part of the ribbons are shown): A, B states are localized 
at the left boundary, with noticeable values of l^^l^ on about 
30 carbon atoms, while l^^l^ of C and D states are completely 
zero in this region. 



states across the width (only part of the ribbons with 
one zigzag edge are plotted). One can observe that the 
wave functions of the A and B states are localized at the 
left boundary, with noticeable values on about 30 carbon 
atoms, whereas the wave functions of the C and D states 
are essentially zero in this region. The contrary occurs 
on the other boundary (not shown in the figure). This is 
topologically distinct from the helical edge states of the 
quantum spin Hall effect, where opposite spins propagate 
in opposite directions along the same boundary. 

The emergence of chiral edge states in the bulk gap is 
intimately related to the topological property of the bulk 
Bloch states in the valence bands. This is characterized 
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by the quantized charge Hall conductance: a.^^ 

where C is an integer known as the Chern number 16|, [l 

and can be calculated from: 
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where f2„ is the momentum-space Berry curvature for the 
n-th band [11, [H, [25,] 



«„(fe) 



-E 






(3) 





FIG. 4: (Color online) (a)-(b) Top view of the spin configura- 
tions of the higher valence band in the kx-ky plane around K 
and K' points. Blue/red regimes represent the spin-up/down 
states. The arrows represent the spin direction at each point. 
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FIG. 3: (Color online) (a) Berry curvature distribution Q (in 
units of e^//i) of the valence bands in the momentum space. 
The first Brillouin zone is outlined by the dashed lines, and 
two inequivalent valleys are labeled as K and K' . (b) The 
profile of Berry curvature distribution along ky — ^. The 
parameters used here are tso — 0.1 and A = 0.18. 



The summation is over all occupied bands below the 
bulk gap, ujn = En/h^ and v^i^y^ is the velocity opera- 
tor. The absolute value of C corresponds to the num- 
ber of gapless chiral edge states along an edge of the 2D 
system [26]. In Fig 131 (a), we plot the Berry curvature 
distribution Q for the valence bands in the momentum 
space. Panel (b) shows the profile of the Berry curva- 
ture along the ky = intersection. We observe that 
the Berry curvature is peaked at the corners of the first 
Brillouin zone. It has the same sign at the inequivalent 
K and K^ points, because the honey-comb lattice pre- 
serves the two dimensional inversion symmetry for which 
Q{k) = Q{—k). This is in contrast with the valley Hall 
effect, found in graphene system with AB sublattice sym- 
metry breaking, where the Berry curvature at K and K' 
points has opposite signs [9]. In the semiclassical picture, 
the curvature field acts on the carriers as a magnetic 
field in the momentum space, and deflects the carriers 
in transverse direction under a driving force. We obtain 
the Chern number from the curvature integration and 
find that C = 2, in agreement with the number of chiral 
edge states from tight-binding calculations. 

There is a simple explanation of the Chern number 
value C = 2. As shown in Fig(TJ(d), there are two valence 
bands: the lower one is almost completely spin-down, 
while the higher one involves a strong variation of the 
spin direction in the momentum space. By plotting the 
spin configuration in kx-ky plane (see FigH]), we find that 
the spin textures around K and K' form two separate 
Skyrmions of the same sign ^] : spins point down in the 
central region (shown in red color), but point up in the 
outside region(shown in blue color). Since each Skyrmion 
contributes to a Chern number C = 1, Skyrmions at K 



and K' points give rise to a total Chern number C = 2 
in our model. 

So far, we have predicted the occurrence of the quan- 
tum anomalous Hall effect for a model Hamiltonian with 
Rashba SOC and exchange field uniformly distributed on 
the honey-comb lattice. In the following, we consider an- 
other model in a similar spirit but with a non-uniform 
distribution: a graphene sheet with magnetized transi- 
tion metal atoms adsorbed on one side (i.e., top). Here 
we study a concrete case of Fe adsorbed graphene us- 
ing first principles method, and results with other kinds 
of adsorbates will be published in a separate paper [28| . 
In our calculations, all atoms are allowed to relax along 
the normal direction of the graphene sheet, but we have 
kept the experimental value of the lattice constant for 
graphene in the x-y plane. Fig IS (a) shows a unit cell 
of a superlattice, with one Fe atom adsorbed on top of 
each 4x4 supercell of the graphene sheet. The calcu- 
lations were performed using the projected-augmented- 
wave (PAW) method [29] as implemented in the VASP 
package [30|, and the GGA exchange correlation poten- 
tial [3l[ was used. 

We find that the most stable adsorption positions are 
on top of the hollow centers (with an equilibrium sep- 
aration of 1.56 A), which is consistent with the previ- 
ous study without SOC [H, |33|. Fig El (b) shows the 
bulk band structure with SOC of the iron taking into ac- 
count. The magnetic moment of iron is found to be about 
2/iB in our calculations. This is consistent with the fact 
that there are totally eight electrons in the 3-d orbitals 
of the Fe atom after two 4-5 electrons being transferred 
there [32|. Therefore, no electron transfer occurs to the 
carbon atoms, thus the Fermi level is exactly located 
in the gap. Figsl5l(c) and (d) zoom in the band struc- 
ture around K and K' points circled in panel (b). One 
can observe that a bulk gap around 5.5 meV is opened 
at both K and K' points. After obtaining the Bloch 
functions from the self-consistent potentials, the corre- 
sponding Berry curvature Q{k) near K and K' points 
are calculated from the Kubo-formula [l6|, l25|: ft{k) = 
— ^^fn^n{k)^ where fn is the equilibrium Fermi-Dirac 
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FIG. 5: (Color online) (a) Typical structure geometry of a 
Fe atom on top of the hollow position of a 4 x 4 supercell of 
graphene. (b) Bulk band structure with SOC included along 
the high symmetry lines. The Fermi level is exactly located 
in the gap. (c)-(d) Zooming in of the band structure at K 
and K' points. A gap ^-^ 5.5 meV circled by the green dotted 
curves is opened, (e)-(f) Berry curvature distribution Q of 
the valence bands near K and K^ points. Note: the Berry 
curvature at K and K' shares the same sign thus giving a 
non-zero Chern number C. 



distribution function. Figsl5l(e) and (f) plot the Berry 
curvature distribution Q summing over all the valence 
bands. Two spikes near K and K^ are clearly seen, and a 
Chern number of C = 2 for the valence bands is obtained 
by integrating over the whole Brillouin zone. This agrees 
with the result of our tight-binding model, despite the 
fact that the exchange and spin-orbit effects produced 
by the Fe atoms are clearly not uniform over the lattices. 

In summary, we find that a non-trivial bulk gap in 
graphene can be produced in the presence of both Rashba 
SOC and exchange field, where we predict a quantum 
anomalous Hall conductance of ayx quantized as 2e^ /h. 
This is followed up with a more concrete example of 
graphene sheet with Fe atoms adsorbed on top. Our 
first principles calculations show a bulk gap as large as 
~ 5.5 meV can be opened at the Dirac points , producing 
the same topological effect. 

Note added: After the completion of the bulk of this 
work (see, the brief announcement in Ref. [34]), we no- 
ticed a similar work [35| demonstrating the quantum 
anomalous Hall effect in the bilayer graphene. 
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